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Abstract 

We consider multiple-input multiple-output (MIMO) radar systems with widely-spaced an- 
tennas. Such antenna configuration facilitates capturing the inherent diversity gain due to inde- 
pendent signal dispersion by the target scatterers. We consider a new MIMO radar framework 
for detecting a target that lies in an unknown location. This is in contrast with conventional 
MIMO radars which break the space into small cells and aim at detecting the presence of a 
target in a specified cell. We treat this problem through offering a novel composite hypothesis 
testing framework for target detection when (i) one or more parameters of the target are un- 
known and we are interested in estimating them, and (ii) only a finite number of observations are 
available. The test offered optimizes a metric which accounts for both detection and estimation 
accuracies. In this paper as the parameter of interest we focus on the vector of time-delays that 
the waveforms undergo from being emitted by the transmit antennas until being observed by 
the receive antennas. The analytical and empirical results establish that for the proposed joint 
target detection and time-delay estimation framework, MIMO radars exhibit significant gains 
over phased-array radars for extended targets which consist of multiple independent scatterers. 
For point targets modeled as single scatterers, however, the detection/estimation accuracies of 
MIMO and phased-array radars for this specific setup (joint target detection and time-delay 
estimation) are comparable. 



1 Introduction 

Inspired by the diversity gains yielded by multiple-input multiple-output (MIMO) communication 

systems for alleviating wireless fading impediments, the concept of MIMO radar systems has been 
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first introduced in [1] and further developed in [2-8] . The underlying idea of MIMO radar for offering 
diversity gains is to illuminate uncorrelated waveforms at the target such that their corresponding 
reflections remain uncorrelated and hence offer independent observations of the target. It has been 
demonstrated that such diversities in observing the target bring about performance improvement 
by enhancing the accuracy in detecting the target and estimating its associated parameters, e.g., 
range, speed, and direction. 

MIMO radar systems can be exploited for capturing radar cross section (RCS) diversity via 
deploying widely-separated antennas [4] , or for establishing more degrees of freedom by configuring 
co- located antennas [9] . While the former configuration improves detection and estimation qualities 
and supports high resolution target localization [5], the latter one enhances the power of probing 
in the vicinity of target location [6, 7, 10]. 

In this paper we consider a widely-spaced antenna configuration and treat two problems. First, 
we analyze the problem of target detection when some radar parameters are unknown and needed 
to be estimated. We offer a framework for joint target detection and parameter estimation (with 
optimality properties provided in Section [3]) when the receive antennas can acquire only a finite 
number of observations. While the proposed framework can be exploited for detecting the target 
in conjunction with estimating any parameter of interest, we consider the problem of detecting 
a target that lies in an unknown location. In our formulation, the uncertainty about the target 
location is accounted for through the time-delays that the transmitted waveforms undergo before 
reaching the receive antennas. We formulate this problem as a composite hypothesis test which is 
shown to be solved optimally via the widely-known generalized likelihood ratio test (GLRT). Note 
that the existing optimality results of the GLRT hold only asymptotically for an infinite number 
of observations under certain assumptions [11, Sec. 5.6]. 

As the second problem, we analyze the diversity gain of the proposed detector, which is defined 
as the decay rate of the mis-detection probability with the increasing signal-to- noise ratio (SNR) 
for a fixed false-alarm probability [4]. We analytically quantify the diversity gain as a function 
of the number of transmit and receive antennas. This notion of diversity gain for MIMO radars 
has been first examined for MIMO radars in [3] and [4] for the MIMO radar model introduced 
in [3], which considers detecting the presence of a target at a given location. It is noteworthy that 
another notion of diversity gain defined as the detector's divergence in high SNR regimes has also 
been analyzed in [3]. 

We treat the two aforementioned problems for two target models: point targets which are 
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modeled as single scatterers [4,5]; and extended targets, which are comprised of many isotropic and 
independent scatterers. The summary of the contributions and results of this work is as follows. 

We propose an optimality measure which is shaped by target detection performance, parameter 
estimation accuracy and false-alarm probability. The corresponding optimal composite hypothesis 
test that satisfies some optimality criteria is introduced and deployed for detecting a target in 
an unknown location. Next, for an Nt x N r MIMO radar system we characterize the maximum 
likelihood (ML) estimate of the time-delay vector which consists of NtN r components corresponding 
to different pairs of transmit-receive antennas. The optimal detector corresponding to such ML 
estimate is also derived. This detector differs from the existing ones in the literature in the sense 
that it aims at detecting a target in an unknown location, whereas the existing MIMO radar 
detectors break the space into small cells and detect the presence of the target in a specific cell. 
Finally, we demonstrate that for the proposed detector, in an Nt x N r MIMO radar system with 
widely-spaced antennas and with an upper-limit constraint on the false-alarm probability, the mis- 
detection probability decays as SNR - ^^'' for extended targets and decays as SNR -1 for point 
targets. Such mis-detection probability in phased-array radars for both extended and point targets 
is shown to decay as SNR -1 . The same observations for extended targets and phased-array radars 
have also been made in [3] and [4], albeit for a different MIMO radar model and different detectors. 

The remainder of the paper is organized as follows. Section [2] provides the MIMO radar system 
model. The statement of the composite detection problem and the definition of optimality are 
provided in Section O We formulate and analyze the joint estimation and detection problem for 
MIMO radars for extended and point targets in Sections U] and respectively, and for phased-array 
radars in Section The analysis on the diversity gain are presented in Section [71 Simulation results 
are illustrated in Section [S] and Section provides the concluding remarks. 

2 MIMO Radar 

We consider a MIMO radar system comprising of Nt and N r transmit and receive antennas, re- 
spectively, and adopt the classical Swerling case I model [12] extended for multiple-antenna sys- 
tems [3,4]. According to this model, a target consists of one or more small scatterers exhibiting 
random, independent and isotropic scintillation. 

We define P as the number of the target's scatterers and denote the locations of these scatter- 
ers in the Cartesian coordination by X p for p = 1, . . . ,P. Also, we denote the reflectivity factor 
of the p th scatterer by ( p and assume that {£ p }^ =1 are identically and independently distributed 
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Figure 1: A MIMO radar system. 



as zero- mean complex random variables with variance p. The target and reflectivity factors are 
assumed to remain constant during a finite number of observations denoted by K and change to 
independent states afterwards. Motivated by capturing the inherent diversity provided by inde- 
pendent scatterers, the antennas are spaced widely enough (such that they satisfy the conditions 
in [3, Sec. II. A]) to ensure that the received antennas capture uncorrelated reflections from the 
target. As illustrated in Fig. [TJ we assume that the transmit antennas are located at X^, for 
m = 1, . . . , N t , and the receive antennas are located at for n = 1, . . . , N r . The transmit an- 
tennas emit Nf narrowband waveforms of duration T with baseband equivalent models given by 
s m (t) for m = 1, . . . ,Nt, where E is the total transmitted energy and J T \s m (t)\ 2 dt = 1. In 
contrast to the conventional phased-array radars which deploy waveforms that are identical upto a 
scaling factor [13], in MIMO radar systems these waveforms are designed such that they facilitate 
acquiring independent observations of each scatterer and often are assumed to be orthonormal [4], 
i.e., 

Vm,n, / s m {t)s* n {t) dt = 5(m — n), (1) 

JT 

where * denotes complex conjugate and <5(-) is the Dirac's delta function. As depicted in Fig. Q3 
the waveform illuminated by the m th transmit antenna to the p th scatterer and received by the 
n th receive antenna passes through an end-to-end distance which we denote by dm,n and undergoes 
a time-delay which we denote by Tm.n = dm,n/c, where c is the speed of light. By defining (3 as 
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the path-loss exponent and superimposing the effects of all scatterers, the base-band equivalent of 
the signal received by the n th receive antenna due to the waveform s m (t) (transmitted by the m th 
transmit antenna) is given by 

p 

r m ,n(t) = \l^rYMli-] s m (t-<„)^" / ' T ™'" +*.>>,»(*)• (2) 



p=l ' 



ation effects of path-losses by including the terms ( -J— ) = ( Tm,n c ) . The exponential term 



Note that this model differs from those of [3] and [4] in the sense that we have added the attenu- 

,/3 / „ \-/3 

exp(— j2?r f c Tm,n) in ([2]) represents the effect of propagation phase shift, where f c is the carrier 
frequency, and z m ^ n {t) ~ CAf(0, jj-), denotes the additive white Gaussian noise. 

We define Xq as the location of the gravity center of the target and denote its associated time- 
delays and distances by {T m n } and {d m ,n}> respectively. We assume that the distances {d m ,n\ are 
considerably larger than the dimensions of the object such that we can replace the distances and 
the time-delays associated with the scatterer X p with those corresponding to the gravity center of 
the target Xq, i.e., eC,n = d m ^ n and r^ n = r m) „,Vp. Therefore 

s m (t- r^J = s m (t- T m>n ), for m = 1, ...,N t , n = 1, . . . , N r , and p = 1, . . . , P. (3) 

Using d2J)- (J3j) , the received signal at the n th antenna, which is a superposition of all emitted wave- 
forms, is given by 

r »W = \ W EE(r) C P s m (t - r m , n ) e~^<« + z n (t) 

m=l p=l ' 

f E- - W)E e~^f^ + z n (t), (4) 




A 

= «T7 



where £ n (i) = Ylm=i z m,n(t) ~ CAA(0, 1). Furthermore, we assume that the waveforms are narrow- 
band. Based on the narrow-band assumption we get [4] 

Vr, s m (t) = e j2nfcT s m {t - r), for m = 1, . . . , N t , and n = 1, . . . , N r , (5) 

which in conjunction with the orthonormality assumption ([1]) implies that 

V T mik , r njfe , / s m (t - T mik )s* n (t - T„ )fc ) dt = 5{m - n). (6) 
Jt 

We also define the time-delay vector r = [t^i, . . . ,Tjv t ,jv r ]- Based on the model given in ([2]) and 
noting that the noise-terms {z n (t)} are unit- variance, the transmission signal-to-noise ratio, denoted 
by SNR, is given by SNR = f . 



3 Problem Statement 



We assume that the receive antennas sample the received signal at the rate of jr samples per 
second. By defining r n [k] = r n (kT s ), z n [k] = z n (kT s ) and s m [A;;r] = s m (kT s — r), the discrete-time 
baseband equivalent of the received signal when a target is present is given by 

r n[k] = J— ^2 -0— h m ,nS m [k;T min ] + z n [k], for k = l,...,K. (7) 

' * m =l Tfn > n 

We also assume that the sampling rate is high enough to ensure that the discrete-time signals 
s m [k;T mt n] remain orthogonal for arbitrary delays T m j,T nt i, i.e., £) fc s m [k; T m ,j]s* [k; T n j] = ^5(m- 
n). Let us define r[k] = [n[k], . . . ,r Nr [k]] T for k = 1, . . . , K and R = [r[l] T , . . . , r[K] T ] T . 

The target detection framework that we propose is different from the conventional target detec- 
tion with MIMO radars. In conventional MIMO radars, e.g., [3] and references therein, the space 
is broken into small cells and the radar detects the presence of the target in a given cell. In such 
radar models, the location to be probed is given, which means that the corresponding time-delay 
vector r is also given. Therefore, the detection problem can be cast as testing whether a target 
exists at a given set of time-delays (c.f. [3]). In contrast, our objective is to detect the presence of 
a target that lies in an unknown location and thereof has unknown corresponding time-delays. For 
this purpose we treat the time-delay vector r as the unknown parameter to be estimated and based 
on that we formulate the target detection problem. In other words we exploit r as an intermediate 
variable which we estimate for being able to carry out target detection. Therefore, in our frame- 
work the main objective of estimating the time-delays is to facilitate performing detection. As a 
minor application, such time-delays can also be exploited for estimating the location of the target. 
Such target localization, however, is not optimal and for achieving optimal target localization one 
needs to treat the target location as the unknown parameter of the interest and deploy the same 
framework for carrying out joint target detection and target localization (instead of time-delay 
estimation) . 

Let fo(R) denote the probability density function (pdf) of the received signal when a target 
is not present. When a target is present, the pdf of the received signal depends on an unknown 
parameter x (in this paper x is the vector of time-delays) and is denoted by fi(R \ x). Therefore, 
by defining the estimate of a; by a;, the detection part of the problem can be cast as 

H : R~f (R), 
Hi: R~h(R\x). 
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Our objective is to jointly characterize the estimator x and find the test for deciding between 
the hypotheses Ho and 7i\. For this purpose we define and optimize a measure which combines 
estimation and detection accuracies. The underlying idea for performing such optimization is to 
deploy the notions of randomized tests and randomized estimators. For any given observation 
vector R £ C'^' we assign complementary probability terms 8q(R) and Si(R) = 1 — d~o(R) to the 
hypotheses Hq and Hi, respectively. In a randomized test, Si(R) for i = 0, 1 is the probability 
that hypothesis Hi is selected. Clearly deterministic tests are special cases of randomized tests 
where {Si(R)} take only the values or 1. Whenever the decision is in favor of Hi, we also 
have to estimate the unknown parameter x. In deterministic approaches, the estimate of a: is a 
deterministic function of the observed data. On the contrary, randomized estimators, instead of 
generating a deterministic value x, generate a random variable x according to a distribution which 
is a function of the observed data. Deterministic estimators can be classified as a special case of 
randomized estimators. 

We strive to identify the decision rules 5q(R) and 5\(R) such that the error in estimating x, 
conditioned on the hypothesis Hi being true, is minimized. This minimization is also subject to 
an upper bound constraint on the false alarm probability Pr(H\ \ Ho)- When the hypothesis Hi 
is true, there exist two types of errors in estimating x. First, the test result can be correctly in 
favor of Hi, but x is estimated erroneously, where we denote the cost function associated with such 
estimation error by C(x,x). Secondly, the test output might be wrongly in favor of Ho, which 
leads to mis-estimating x and results in an error cost which we denote by Co (a;). Note that the 
cost function C(x,x) can be selected arbitrarily. In this paper we consider the mean-square error 
(MSE) cost given by C(x,x) = \\x - x\\ 2 and Co(x) = \\x\\ 2 . 

Subject to a maximum tolerable level of false-alarm probability Pf a , the estimation strategy 
which minimizes the average cost over all possible randomized estimators is given by 

I min {Si(R)} K R,x,x [C(x, x) | Hi] + Er )X [C (x) I Hi] 

r{R) - < (yj 

^ s.t. Pr (Hi | Ho) < P fa . 

The optimal hypothesis test (formalized by {5i(R)}) and the optimal estimation strategy yielded 
by solving the problem V(R) for the choices of C(x, x) = \\x — x\\ 2 and Cq(x) = \\x\\ 2 are described 
in the following theorem. 

Theorem 1 (Moustakides [14]) For a finite cardinality vector R G Cl H l and an unknown vector 
parameter x £ A, solving V(R) provides the optimal estimation strategy x and the subsequent 
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optimal detection rule for deciding between TIq and H\ are given by 

fi(R I x) w i 

x = argmax7r(a;)/i(i? | x), and — — ^ A, (10) 

where ir(x) is the prior distribution of x and the estimation cost functions are C(x,x) = \\x — x\\ 2 
and Co(x) = \\x\\ 2 . 

The above theorem essentially establishes the GLRT as the optimal estimation/detection strat- 
egy with properties P1-P4 enumerated in the sequel. Several other asymptotic optimality results 
are known for the GLRT which are all based on having an infinite number of observation (c.f. [14] 
and [15, Chapter 22]). x given in (jlOh is the maximum a posteriori (MAP) estimate of x. As we 
do not have any prior information about the parameter to be estimated (the vector of time-delays 
r), throughout the analysis we assume that ir(x) has a uniform distribution. Hence, the MAP 
estimate of x becomes its ML estimate. The estimator and detector provided by Theorem [1] have 
the following properties. 

PI) The ML estimator x minimizes the average minimum- mean square estimation error of x. 

P2) The false-alarm probability of the target detector is kept below a certain level. 

P3) For the given set of ML estimates x, the target detector is Bayesian-optimal, i.e., the Bayes 
risk is minimized [16, Sec. II. B]. 

P4) The test requires only a finite number of samples, i.e., \R\ < oo. 

In order to apply the theorem above to our joint estimation/detection problem we set the vector 
of time-delays as the parameter to be estimated, i.e., x = r. For deploying the optimal estimator 
and detector provided in (|10p we also need to characterize the domain of x = t which in Theorem 
[I] is denoted by A. Note that the vector r cannot be any arbitrary vector in the space M^ tNr . 
This is due to the fact that for any given placement of antennas, there exists a correlation between 
the elements of r and there are some equations relating these time-delays. For the given antenna 
placements {^m}m=i anc ^ {^n\n=i ^ us define X = {A"^,A^}. Also for the given antenna 
placement X define A(X ) C W^_ tNr as the set which contains the vectors that can be valid choices 
for r. 

Precise characterization of A(X) is crucial for attaining the optimal estimation and detection 
performance. Any inaccuracy in such characterization degrades the performance and provides only 
sub-optimal estimators and detectors. In the sequel we provide necessary conditions for a given r to 
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be a valid time-delay vector. Such necessary conditions describe a set of time-delay vectors which is 
a superset of A{X). We also conjecture that these conditions are sufficient which consequently can 
characterize A{X) precisely. If these conditions are not sufficient, then the resulting approximation 
of A(X) will only provide a sub-optimal solution. 

For establishing the necessary conditions let us break down the end-to-end time-delay T m ^ n as 
T~m,n = tm + t' n where t m is the time required for the waveform to travel from the m th transmit 
antenna to the target and t' n is the time required for the reflected waveform to travel from the target 
to the n th receive antenna. If a given r is a valid time-delay vector then Vm,i € {1, . . . ,Nt} the 
three terms t m , ti, and - — ^ — — should constitute the lengths of a triangle with its vertices given 
by the locations of the m th transmit antenna, the i th transmit antenna, and the target (||A^ — X\\\ 
is the distance between the m th and i th transmit antennas). The triangular inequalities impose 
that r is a valid time-delay vector only if 

\t m ~k\ < ll X ™~ X *H for m,i=l,...,N t . 

c 

Similarly, by considering the triangles corresponding to the receive antennas we get the necessary 
conditions 



14-4-1 < " _ 3 " for n,j = l,...,N r . 



c 

By recalling that r m>n = t m + t' n V m,n,i,j, we find that \t m — ti\ = \r m ^ n — r^ n \ and \t' n — t'^ 
\T m ,n — T m,j\- Therefore, for the subspace A(X) we have 



f ILX* ||X f r t -X r || ] 
A{X) <Z It : \T m)n - n : n\ < — ^— — l — , and | T m ,n - T m ,j | < — — V m,i, j,n *> . (11) 

We conjecture that any vector r that satisfies the conditions above is a valid time-delay vector 
and therefor the subset operator C can be replaced by equality. We underline that, however, in 
the case that these conditions are not sufficient, performing optimization on the above superset of 
A{X) will provide only sub- optimal estimators and detectors. 



4 Joint Detection and Estimation for Extended Targets 

In this section we consider the targets that are extended enough to be modeled as a group of 
isotropic and independent scatterers, i.e., P > 1 in ([2]). Our objective is to optimally detect a 
target, based on the model in ([8]), and simultaneously estimate the vector of time-delays r. At 
the same time we assume that only a finite number of observations are available. According to 
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Theorem Q] this problem can be solved via the GLRT given in ()10p for x = r and for .A = A(X) 
given in (jlip . 

We start by deriving the ML estimate of t and then provide the detection-related analysis. 
We also develop a technique in Section [8] based on the Newton- Raphson method for localizing the 
target by using the estimated time-delay vector. Exploiting this simple target localization method 
is not the optimal strategy and is merely intended to provide some numerical evaluations to show 
the improvement gained by deploying MIMO radars. 

4.1 Time Delay Estimation 

The hypothesis given in ([8]) can be cast as 
Ti-o ■ r n [k) = z n [k], 

fW N i ( 12 ) 

T~L\ : T n [k] = \ Jq~ C zZm=l ~~B ^m,n Sm^i^m,!)] ~\~ Z n [k]. 

Define h = [hn, /ii,2, • • • , hN t ,N r ] T i where h mn , as defined in accounts for the effects of the 
position and reflectiveness of the scatterers corresponding to the m th transmit and the n th receive 
antennas. The antennas are widely separated and the reflectivity factors {Cp} are complex and 
independently distributed with zero mean and variance -p with P>1. Therefore, by using the 
central limit theorem, {h mtn } are i.i.d. and distributed as CM (0,1) [3]. 

Based on the model given in 0, when a target is present, the vector of the received signals 
R depends on the time-delays, which we are interested in estimating, as well as on the unknown 
random vector h. Therefore, when a target is present, the pdf of R for any given r and h is 
fi(R | T,h). In order to obtain the ML estimate of r one approach is to estimate it through solving 
arg max Tg _4(^) fi{R \ t), which requires recovering f\(R \ r) from f±(R \ T,h) by averaging over 
all realizations of h. Alternatively, we can estimate it jointly with h when doing so is deemed to 
be beneficial. The latter approach is more beneficial when an accurate estimate of h is available, 
e.g., in the high SNR regimes, while averaging leads to a better performance when the estimate is 
very inaccurate, e.g., in the low SNR regimes. The ML estimate of r is provided in the following 
theorem for both scenarios. 

Theorem 2 The ML estimate of r for extended targets and a given antenna placement X 
1. through MAP estimation of h is given by 



N t N r 



t = arg max > > 1 T 

MV\ ^ t-^ 1 



X^fc=l r n[^] s m[k'i T rn,r 
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2. and through averaging over all realizations of h is given by 

2 



Aft N r 

arg max > > 

m=l n=l 



T„ ^ E \^'m,n) 



log 



>cr, 



Proof: See Appendix lAl 



4.2 Target Detection 

Based on the ML estimates of the time-delay vector r provided in Theorem [21 we proceed to find 
the corresponding optimum detectors. We show that both estimates give rise to the same optimal 
detector given in the following lemma. 

Lemma 1 The optimal test for extended targets and for the given estimate f is 

2 



Nt N r 

EE 

m=l n=l 



TsNi 



Hi 
Wo 



(13) 



Proof: Based on the model in (|12p the likelihood fo(R) under the null hypothesis is given by 



fo(R) = (^r^exp^-^llrMII 2 ^. 



(14) 



Finding the estimate -r MAP is carried out by jointly estimating r and h. For this estimate of r, by 
setting x = [r, h] the optimal detector characterized in (jlOh is found as 

m&x x&A ir(x)fi(R | x) _ max TE-4( ^ )ife f x {R \ r,h)f h (h) _ f^R \ T,h)f h (h) Hx 

fo(R) ~ MR) ~ MR) Ho ' 

where r and h are the MAP estimates obtained in Theorem [2j By following the same line of 
argument as in Appendix [S] ([30]) - ([39[) and recalling the definitions of a T and B T given in ([32]) and 
([35]) . respectively, the detector is given by 



log 



h{R\T,h)f h (h) 

MR) 



log 



■ (tt)-^ exp (- £ fc ||r[fc]|| 2 ) exp (a£l*T V) 



(vr) 



N t N r 

EE 

m=l n=l 



2/i 



Wi 
Wo 



where 6* = ^ log A. We further define b m ^ n = Y^k=i r *nW\ s m[k; f mjV ] and l m>n = + (cf„ 
Therefore, the value of the threshold 9 is linked to the probability of false alarm Pf a as 

IK, 



(15) 
(16) 

\2f3 



p f = p EE 



^ m=l n=l 



> 
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Note that under 7io, &m,n is distributed as CJ\f(0, jr) and therefore, |6 mj n| 2 is exponentially dis- 
tributed with parameter T s . Hence, {|6 m ,n.| 2 /^m,n} have exponential distributions with distinct 
parameters T s ■ l^ n . Thus, Yl m Yin * s sum °^ NtN r independent exponential random 

variables with different parameters and its pdf is given by [17] 



j£ ^_ ^""^ ^""^ \ u m,n\ 



\K ' 2 



,m=ln=l 



Aft N, 



T ■ 1 



N t N r _ / rp ,-1 



ex P ( 2s lm,n x ) 
m =l n =l EL^n T s(^ij ~ l m,n) 



If we denote the cumulative distribution function (CDF) of this random variable by Gi(-), then 9 
can be determined as = G^ 1 (1 — Pf a ). By following a similar line of analysis for the case that we 
deploy the estimate r ave , the same detector is found. ■ 
It is noteworthy that our detection scheme has two major differences with that of [3] provided in 
[3, Eq. (24)]. First, the detector in [3, Eq. (24)] tests whether a target exists at a known time-delay 
while we try to detect a target with unknown time-delays. Secondly, the model of [3] embeds the 
effect of the time-delays only as phase shifts and ignores the path-loss effect. By recalling that the 
path-losses also depend on the time-delays we have modified the model of [3] to also capture the 
effects of path-losses. 



5 Joint Detection and Estimation for Point Targets 

In this section we consider the application of MIMO radars for detecting point targets or exposing 
single-scatterers. In this case, the target is modeled by one scatterer located at its gravity center 
X . Note that in ((!]) by setting P = 1 we get h m>n = £ e ~ j27rfcTm - n . Unlike in extended targets, 
the previous section, the distribution of the reflectivity factor £ is not known and thereof it should 
be treated as an unknown quantity to be estimated along with the time-delays. Therefore, the null 
and alternative hypotheses are given by 

{H : r n [k] = z n [k], 

As before, we are interested in solving the optimum test given in (jlOp and we start by determining 
the ML estimates of the reflectivity £ and the time-delay vector r. 

5.1 Time Delay Estimation 

Based on the model given in (|18p and the fact that the distribution of £ is not known a priori, we 
estimate r along with £. The pdf of the received vector R for any given r and £ is fi(R | t,£). 
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By setting x = [r, Q and using the optimal test given in (jlOp the ML estimate of r can be found 
as shown in the following theorem. 



Theorem 3 The ML estimate of r for point targets and antenna placement X is given by 



t = arg max 
t(EA{X) 



m=l Z_m=l 



1 m , n 



sr^N t v^7V r 1 



1 m,n 



Proof: See Appendix [Bl 



5.2 Target Detection 

Based on the ML estimate of t provided in Theorem [3] the optimal detector is characterized by the 
following lemma. 



Lemma 2 The optimal test for point targets and for the given estimate t is 

V t N r -Q-Trf t K 2 

Tm.n 



Hi 
Ho 



m=ln=l ' m >™ fc=l 

Proof: As we are estimating r jointly with by setting x = [r, £] an d applying Theorem [1] the 
optimal detector is given by 

max xeA Tr(x)f 1 (R \ x) = max TeA{x)jC frjR | t,Q = h(R | t,Q * i . 

/o(i«) /o(fl) /o(H) Wo ' 

By following the same line of argument as in Section 14.21 and by using the estimates r and £ given 
in Theorem [3] and (|44p , respectively, the optimum test for a point target is given by 



log 



h(R\r,C 
MR) 



log 



(vr) 



-KN, 



\ 



expC-E.HArimexpf^lCl 2 £^ =1 En=i J?) 
(vr)-^-exp(-E fc l|r-M|| 2 ) 



\ 



^m=i Ln=i ( CT m,n> log A, the test can be cast as 



„-2/3 
72T3- 

' m , n 



Moreover, by setting 



Hi 

Z A. 

Ho 



(19) 



m=ln=l Tm > n k=l 



Hi 
Ho 
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In order to determine the value of the threshold 9, note that e j27r ^ cTm '" 5^fcLi r n[k] s^[k; f m n ] is 
distributed as CM(0, jr) under Ho and X)fcLi r n[k] s^[k; f min ] is independent of Ylk=i r n' [k] s* m , [k; f m f >n i 
for m / m' or n / n'. As a result, under 7i$ we have 

3. ^ ^/.W * / 1 ^ ^ 1 \ 

— ^ Z_^ r n[k] s* m [k;T m , n ] ~ CAA ( 0, — 2^ 2^ ^ 



m=l n=l 



fc=l 



2/3 

m=ln=l 'm,n 



and consequently 



K 



EE 

m=l n=l 



1 m,n 



fc=l 



Exponential 



Z_/m=l Z_m=l 



iVr 1 



f 2/3 
' m,?i 



Therefore, for a given value of Pf a , the threshold level 6 is found as 



N t N, 



s m =l n=l Tm > n J a 



6 Optimal Joint Estimation/Detection for Phased-array Radar 

We will compare the performance of MIMO radars against that of conventional phased-array radar 
systems. A phased-array system utilizes an array of closely-located antennas and therefore, the 
channel coefficients h m ^ n are fully correlated, i.e., h m<n = h. For the case of extended targets h 
is distributed as CJ\f(0, 1) and for point targets its distribution is unknown [3]. Moreover, all the 
emitted waveforms are equal up to a scaling factor, i.e., s m (t) = s m s(t) for m = 1, . . . , Nf such that 
Ylm=i \ sm\ 2 = Nt (for having the total transmitted energy equal to E). Therefore, by using (|7|) 
and the narrow-band assumption, the system model is given by 

FT Nt l 

r n [k) = J w c~Ph s[k;r ltl ] £ — S me J2*/°Ki-T,»,») + Zn[k]: {ot k = 1, . . . ,K. (20) 

V * m=1 Tm,n 

The following lemma is instrumental to deriving the ML estimate of the time-delay vector. 
Lemma 3 For any given set of functions {gi{t)}f = i where gi(t) : M — > C and a,fGE, we have 

N 

= X>(*i)i- 

1=1 

Proof: See Appendix O ■ 



max 
t 2 ,...,t N & 



N 



E eiati ^(*i) 



i=l 
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For the purpose of analyzing the diversity gains as given in Section [7J as well as comparing es- 
timation performance, we provide the optimal detector and estimator for phased-array radars in 
the following lemmas. For further analysis we define 5 n (r) = Y^*_-i -, — ^-^s m e^ 2 ' K ^ c ^ Tl ^~ Tm ^ and 

'ill' J- yCTm^n) 1 

SJ T ) = , 1 ^ s m e j27rf ^ Tl ^- 2T ^l Similar to A(X) which we defined for MIMO radars, 

we also define A'(X') such that it contains the valid choices of the time-delays for a phased-array 
radar with closely- located antenna placement X' . 



Lemma 4 For a given antenna placement X' , the ML estimate of the time-delay vector r and the 
optimal detector in phased-array radars are as follows. 



1. For extended targets through estimating h we have 



~ MAP 

r = arg max 



£n=i^(r)££ = iM^;n 



K 



kY%U\Sn{T)\ 2 + l 



2. For extended targets through averaging over all realizations of h we have 



En=l S n( T ) £fc=l r n[k]s*[k; ri,l] 



T7 Y1n=1 I^( t )| 2 + 



log |^X> n (T)| 2 + l 



T s N t 



n=l 



3. For point targets through estimating £ we have 



t = arg max 
-reA'(X') 



En=l Sn(T) Efc=l r n[k]s* [k] n,l] 



Proof: See Appendix [Dj ■ 
Based on the ML estimates of r provided in Lemma [U the optimal detectors are presented in 
the following lemma. 

Lemma 5 The optimal test for the given estimate r is given as follows. 



1. For extended targets 



2. For point targets 



N r 



K 



J2 S n(r)^2r n [k]s*[k;f 1A ] 



71=1 k=l 



Hi 
Ho 



N r 



K 



71=1 fe=l 



2 Hi 
Ho 
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Proof: For the case that we are estimating r and h jointly, we set x = [t, h] and therefore, the 
optimal detector is given by 



max xeA it (x) MR | x) _ max^^g | r,h)f h (h) _ f^R \ r,h)f h (h) Hi 

MR) ~ MR) ~ MR) no 



% A. 



Proceeding with the same line of argument as in Section 14.21 the optimal test is given by 

h(R\T,h)f h (h) _ \ \ZtlS*n(T)Zk=irn[k} S *[h,f hl }\ 2 } Hx > 

MR) exp l ^£^l^)l 2 + f J ho ' 

or equivalently 

N r K 2 H 

^2S* n (f)Y,rn[k}s*[k;f 1A 

n=l k=l 

where 

e = io g x.(^f:\s:(r)\ 2 + § 

^ s n=l 

Under Ho, the term S*(r) Ylk=i r n[k] s* [k; f\j\ is distributed as CJ\f(0, ) and ^fc=i r nW s *[^5 ^1,1] 

is independent of J2k=i r n'[k] s*[k; f^i] for n 7^ n'. As a result, under we have 



Ho 



N r K / N r \ 

£ S£(f ) £ r„[fc] f liX ] ~ CAT 0, - £ |5 n (r)| 2 

n=l k=l V s n=l / 



and consequently 

N r K 2 



n=l fc=l 



Exponential 



En=l\Sn(T)\ 2 



Therefore, for a given value of Pf a , the threshold level 6 is found by 

1 ^ 1 
9 = -Y,\Sn(r)\ 2 log p-. 

We can similarly find the same detector when the estimate T ave is computed by averaging out the 
effect of h. ■ 



7 Diversity Gain Analysis 

In order to compare the performance of MIMO and phased-array radars quantitatively, we analyze 
the diversity gain, which determines how fast the mis-detection probability decays as the trans- 
mission SNR increases while there is an upper-bound constraint on the false-alarm probability [4]. 
More specifically we are interested in finding 

.. logP md (SNR) 
a = — nm — - — — — — , 
SNR^oo log SNR 
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where P mf ;(SNR) denotes the mis-detection probability at the signal-to-noise ratio SNR. The di- 
versity gain for the MIMO radar model provided in [3] has also been examined in [3,4]. [3] also 
analyzes another notion of diversity, which is the asymptotic value of the detectors's divergence. 
Throughout the analysis in the sequel, we say two functions fi(x) and f2(x) are asymptotically 

denote it by f\(x) = f 2 (x). 



equal when lim^oo = 1 and denote it by f\{x) = f2(x). We also define < and > accordingly. 
We also say two functions fi(x) and f2(x) are exponentially equal when linx^oo j° s £ p) = \ an d 



Remark 1 Asymptotic equality is a sufficient condition for exponential equality. 

We will use the following lemma for analyzing the diversity gain of MIMO and phased-array radars. 

Lemma 6 For any M independent Gaussian random variables Y m ~ M(p- p m , cr 2 ), m = 1, . . . , M, 
where p G M + and p m ~ M(0, c^J, and for any given 7 E M + , in the asymptote of large values of 
p we have 



/ M \ 

Pr [J2 Yl < 7 



\m=l 



P~ M , (21) 



where /x = [pi, . . .,p M }- 



Proof: See Appendix [El ■ 
By using the lemma above, in the following theorem we establish the diversity gain achieved by 
MIMO and phased-array radars for both extended and point targets. -P m d(SNR). 

Theorem 4 The diversity gain achieved by 

1. an N t x N r MIMO radar system for extended targets is N t N r , i.e., P^ d (SNR) = SNR _7Vt7Vr ; 

2. an N t x N r MIMO radar system for point targets is 1, i.e., P^ d (SNR) = SNR -1 ; 

3. an Nt x N r phased array system for extended and point targets is 1, i.e., P_?4(SNR) = SNR -1 . 
Proof: 

1. For extended targets, by recalling the definitions of b m ^ n = J2^=i r nW\ s m[k; T mj „] and l m>n = 
T^Ni + (cT m)n ) 2/3 , from the optimal test in (fT3l) we have 



P^(SNR)=Pr(^^i^< 



1 1 lm,v 
\m=l n=l ' 



Hi . (22) 



17 



We define b 7 mn = 3?(& mi n), b l mn = ^s(b mn ). Therefore, the mis-detection probability given in 
is equivalently given by 

N t N r 



P^(SNR) = Pr J- 

\m=l n=l ' 



< 



Hi . 



(23) 



By further defining l m \ n = min min l m ,n and Z max = max m>n l m>n we can provide the following 
upper and lower bounds on the probability in (|23l) 



/ N t N r 



yn=l n=l 



Wi <^ d (SNR) 



/ iV t N r 



\m=l n=l 

On the other hand, under hypothesis H\ from (112j) we have 



Hi 



(24) 



FT Nt 1 K K 

* m'=l ^~m'n k=l k=l 



which in turn provides 



^-8(m— m') 




Furthermore, by defining h r mn = lSt(h mjV ) and h l mn = ^s(h m:n ) we have /ij^^, fa^n ~ A/"(0, |) 
and 



h r 

-P '"m,n 



1 



(cr m , n )^.^,— I and fe^-JV 




T s ' 2T S 



Now, we apply Lemma [6] by setting M = 2N t N r , {Y m } m = {b r ,U ' } m , n , P = v#, 



{^m}m = a/ 7^(cr mi „) ^ x {h r mn , h l mn }, a 2 = ± and = \ . Thus, we get V7 



/ Nt N r 

Pr E E( W + 

\m=l n=l 



Hi = (V£) 



-2N t N r 



(T s SNR)- NtNr = SNR 



-AT t Af r 



Hence, for the choices of 7 = 6l m \ n and 7 = 6>/ max we find that the lower and upper bounds 
on P^ d (SNR) given in (|24p are both exponentially equal to SNR -ArtAr, \ This results in an 
identical exponential order for P^ d (SNR), i.e., 



P^(SNR) = SNR 



-NtN r 



which establishes the desired result for MIMO radars with extended targets. 
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2. We define 

* - EE 



~y ~] r n [k] s m [k; T m ^\ 



m=ln=l m ' n fc=l 



m=l n=l Tm >™ 



3. JE, fjZKfcfnr K 



m'=l rri'.n fc=l 



Tj^5(m—m') 



EE 



/9 



=1 n =l Tm >™ fe 



Therefore, 6 is distributed as 



\ V JV * m=ln=l IWWJ 1 s m=ln=l T s T m ,nJ 

By following a similar line of argument as for extended targets and defining b r = R(6), 
b l = 3(6) and 

JE, JE, gj27r/ c (f m , n -r m , n ) \ 



a fJ^Sk e i27r/ c (f m , n -r m ,„) 

EE 



1 n=l ' ' / \m=l n=l ' ' 



c 



we get 



% N r 



r„tf(JE c -fi.rYY*^) t an d &~tflJE c -fi.CYY — 



m =l n =l 2T s T m , n 



TE ^ ^ 2T f 2/3 

s m=l n=l zl sTm,n 



and 



P^(SNR) = Pr (|6| 2 < <?J = Pr ((^) 2 + (6*) 2 < 

Now, we apply Lemma[6]by setting M = 2, {Y m } m = {b r , p = v 7 ^, {/i m } m = A / 4-c _/3 x 
{C r ,C}, ^ 2 = ^ an d <*i = E[(C) 2 ], of = E[(<*) 2 ]. Therefore, 



P^(SNR) = (VE)- 2 = (T.SNR)- 1 = SNR 



-i 



3. We provide the proof for extended targets. Similar line of argument and appropriate modifi- 
cations provide the proof for point targets. Based on the optimal test for detecting extended 
targets by phased-array radars given in Lemma [5] we have 

/ N r K 



pPA _ p 



n=l 



k=l 
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2 




< 9 





Under hypothesis Hi from (|20p for extended targets we have 

N r K rjr N r K 

6 = EW)E r »W s *[*i^d = V/ F ^^ (f)5n(T)/l ^ s[fe;T]s * [A:;f] 

n=l fc=l n=l k=l 

s .. ' -> ' 



Af J_ 
= h T s 



N r K 



n=l k=l 



CJV(0,^) 



where /i ~ CAA(0, | ^^=i 5*(t)5„(t)| 2 ). Therefore 




E h N r 



N t T S 'T S J- 

Similarly as before we define b r = U = /i r = and /i* = 9(A), where we have 

h r , U ~ CAA(0, IE "- lS y )S " (T)|2 ) and 

and ^C7vf#.|,^V 

Again we apply Lemma [6] by setting M = 2, {Y\,Y2\ m = {b r ,b}, p = \/~E, {jUi,/i2}m = 
-^(ct)- 13 x {/i r ,^}, a 2 = and cr^ = Thus, we get V 7 

Pr ((6 r ) 2 + (b 1 ) 2 < 7) = (V^)" 2 = SNFT 1 . 

As a result we find the following diversity gain achieved by phased-array radar systems for 
extended targets 

^(SNR^SNR- 1 . 

■ 

Although we are using a MIMO radar model different from [3] and subsequently derive a different 
detector, the results above conform with those of [3] which considers detecting the presence of a 
target at a given location. This result demonstrates that in terms of diversity gain, our proposed 
MIMO radar model is capable of capturing the same diversity gain achieved by the MIMO radar 
model of [3]. 



8 Simulation Results 

8.1 Extended Targets 

In this section we provide simulation results on the performance of the proposed joint estima- 
tion/detection framework. We consider two antenna configurations with Nt = N r = 2 and 



20 
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Figure 2: Average normalized MSE of time-delay estimates versus SNR for extended targets. 

Nt = 4, N r = 8. For the MIMO radar we assume that the transmit and receive antennas are 
located at = (m, 0, 0) for m = 1,...,4 and X r n = (0, n,0) for n = 1,...,8, respectively. 
For the phased-array radar we assume that the transmit antennas are all closely-located around 
X l m = (1, 0, 0) for m = 1, . . . , 4, and the receive antennas are closely-located around X r n = (0, 1, 0) 
for n = 1, . . . , 8. Also we assume that the target to be detected is located at Xq = (20,15,0), 
where all the distances are in kilometer (km). The path loss coefficient is (3 = 2 and the carrier 
frequency is f c = 5 MHz. We assume that the target comprises of P = 10 scatterers and the 
number of signal samples is K = 40. Finally, for the MIMO radar the emitted waveforms are 
s m (t) = exp ( ^ 2 ™ lt ^ (U(t) —U(t — T)), where U(t) is the unit step function and T denotes the 
duration of the waveform and the sampling rate at the receiver is T s = jq. For the phased-array 
radar all the emitted waveforms are equal to Si(t). 

We first consider the performance of parameter estimation. Fig. [2] depicts the average normal- 
ized mean-square errors (MSE), i.e., N ^ N ^2 m Y2n \ Tm " Tm '" ^ or phased-array and MIMO radars 
as a function of received SNR. It is observed that the MIMO radar outperforms the phased-array in 
all SNR regimes and in particular, by a large margin in the low SNR regime, which is of more interest 
in radar applications. Moreover, it is seen that in MIMO radars, the MAP estimator f- MAP per- 
forms better than the estimator r ave in the high SNR regime, while the estimator r ave outperforms 
the MAP estimator r MAP in the low SNR regime, as expected and discussed in Section 14.11 
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Figure 3: Average normalized MSE of location estimates versus SNR for extended targets. 
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Figure 4: Probability of mis-detection versus SNR for a false alarm Pf a = 10 6 for extended targets. 

It should be noted that comparing the accuracies of the time-delay estimators are not as in- 
formative as comparing the accuracies of the detectors. This is mainly because the time-delay 
estimates have different roles in MIMO and phased-array radars. Nevertheless, we perform such 



22 



— b — MIMO, N=N =2 

t r 

o MIMO, N=N r =4 

_□_ PA, N=N=2 

t r 

_o_ PA, N=N=4 

t r 

°i 1 0^° KT 8 KT 6 10" 4 10" 2 10° 

Probability of false alarm 

Figure 5: Probability of target detection versus probability of false alarm for SNR = dB for 
extended targets. 

comparisons as the detectors accuracies in our formulations depend on the estimators accuracies, 
which justifies comparing the time-delay estimator in MIMO and phased-array radars. 

We next illustrate the localization performance of MIMO radar in Fig. [3l Once the vector r is 
available, it is then possible to further estimate the location Xq of the gravity center of the target. 
The localization problem can be formulated in the form of the following nonlinear equations 

f = <j>(X ) + w, (25) 

where u is the noise term, and 4>(Xq) = [<pi t i(Xo), . . . , <f>Nt,N r (Xo)] is the vector of functions defined 
as 

(pm,n( X o) = ^{\\X -X t J\ + \\X -Xr\\),m = l,...,N t , n = l,...,N r , (26) 

where X^ and X^ are the positions of the m th and n th transmit and receive antennas respectively. 
We can solve for Xq from ([25]) iteratively as follows. Denote X^ as the solution at the i th iteration. 
By linearizing <^>(Xq* +1 ' ) ) with respect to Xq we get 

<t>(xt 1] ) « + H (X^) (Xt 1] ~ 4h ( 2 7) 

where Hq(Xq ) is the Jacobian of the function </>(.) evaluated at Xq . Then, the least-square 
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estimate of the position Xq based on (|25|) and (j27j) is given by 



v {i+l) 

Xq = argmm 



y (X®)-H (X®)X y{X*>) - H (X*>) X , (28) 



T 



(*)\ u~„ / v(i)^ 



where y(X^) = r — ( cj)(X^) — HqXq^ \ and therefore Xq* +1 ^ has to satisfy the normal equation 



H (X®) T H (XW)\ X^ = HoWfvW)- (29) 

An initial estimate X^ can be obtained, e.g., using the method in [18]. It is seen that the 4x8 
MIMO radar performs considerably better than the 2x2 MIMO radar, which is due to the fact that 
32 time-delays provide much more information about the position of the target than 4 time-delays 
do. 

Finally we consider the detection performance. In Fig. H]the probability of mis-detection versus 
SNR is illustrated. The tests are designed such that the probability of false alarm is Pj a < 10 -6 . As 
analyzed in Section [7] and observed in this figure, the slope of the mis-detection probability of the 
phased-array radar is 1 decade per 10 dB, whereas that of the MIMO radar is NtN r times steeper. 
Figure [5] shows the receiver operating curve (ROC) for the MIMO and the phased-array systems, 
for SNR = dB. It is seen that the MIMO radar significantly outperforms the phased-array radar 
over a wide range of false alarm values. 

8.2 Point Targets 

We use a similar system setup as in Section [8.11 Fig. [6] illustrates the normalized average MSE in 
estimating the time-delays versus SNR. Unlike for extended targets, it is observed that for point 
targets, conventional phased-array and MIMO radars exhibit similar target detection and time- 
delay estimation accuracies. Therefore, when considering joint target detection and time-delay 
estimation, deploying MIMO radars are not much advantageous for point targets. This conclusion 
is nevertheless limited to the specific problem analyzed in this paper and MIMO radars can be 
potentially advantageous in other scenarios like those discussed in [4] and references therein. This 
is due to the fact that point targets lack independent scattering section and thereof provide no 
diversity gain. 

In Fig. [7] and Fig. [5] the probability of mis-detection versus SNR and the ROC are plotted, 
respectively. For Fig. [TJthe tests are designed such that the probability of false alarm is Pf a < 10~ 6 
and for Fig. [5] we have set SNR=0 dB. It is seen from Fig. [7J that both the phase-array and the 
MIMO radars exhibit a diversity gain of 1, which verifies Theorem HJ 
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Figure 6: Average normalized MSE of time-delay estimates versus SNR for point targets. 
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Figure 7: Probability of mis-detection versus SNR for point targets and Pf a 
targets. 
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Figure 8: Probability of target detection versus probability of false alarm for point targets and 
SNR=0 dB for point targets. 

9 Conclusions 

In this paper we have offered a framework for jointly detecting the presence of a target and es- 
timating some unknown target parameters. We have provided a composite hypothesis test that 
satisfies some optimality criteria for both target detection and parameter estimation. By using 
this framework we have proposed a new MIMO radar model for detecting the presence of a target 
in an unknown location and formulated a composite hypothesis test for solving this problem. In 
this problem, the unknown parameter of interest to be jointly estimated along with target detec- 
tion, is the vector of time-delays that a transmitted waveform experiences from being emitted by 
the transmit antennas until being received by the receive antennas. For the proposed detection 
scheme, we also have analyzed the decay rate of the probability of mis-detection with increasing 
SNR while the false-alarm probability is kept below a certain level (diversity gain). Our simulation 
results demonstrate that for the specific problem of joint target detection and time-delay estima- 
tion, deploying MIMO radars is advantageous for extended targets. On the other hand, for point 
targets, the accuracies of target detection and time-delay estimation by conventional phased-array 
are comparable with those yielded by MIMO radars. 
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A Proof of Theorem [2] 



The pdf of the vector of received signals R for any given r and h is given by 



N r K 



f 1 (R\r,h) = (n)- KN >- e xp[ -£E 



n=l k=l 



He Nt 1 

V t -1 Tin r> 



m=i ' m > n 



(30) 



=Dl(x,/i) 

By taking into account the orthogonality assumption i s m \k\ T mM 
can further simplify the term D\{r, h) defined in (|30p as follows. 



Af r A 



Dl(T,fc)=££ 



n=l fc=l 
A 



r^- 1 

V I 1 Trr) rj 



AT r N t N t 



^2 ll r W II 2 + c ^ 5^ "9 /3 h m \n s m[fe; T m,n] s m / [^' T m / ,n] 



fe=l 



T /3 J3 
n=l m=lm'=l m > n m' ,n 



fe=l 



7^- 6(m—m') 



Fe Nr Nt l K 

V Jj~ ° 13 £ £ ^ m ' n — fl £ r nM s m[^5 T m,n] 

V * n=lm=l Tm .« fc=l 

f-g- N t 1 A 

\ C ^ y / ^Ti>" ~ S ^1 r w [fc] s m[^i T m,n] 

V £ i 1 Trr? r) 1 



n=l m=l 



fc=l 



fe=l 

N r N t 



T s N t — ' —i ; 

n=l m=l 



n=l m=l * Tm > rl fc=l 



a m,n( T ™.™) 

A' 



^ ^ c -/3 

^m,n \l ~J3 y> r nW s *m[^ T m,n]- 

n=l m=l 



(31) 



Q>m,7i {j~m,n ) 

For any given time-delay vector t, define 
, a E 



T s N t 



c ^diagfa ™,r, 2 , -,t n % J and a x = [01,1(1-1,1), ai, 2 (ri )2 ), . . . , a NuNr (T NuNr )} J 



Then, Di(t, h) can be represented in the matrix form as 



A 



Dx{T,h) = ^2\\r[k}\\ 2 + h H A T h-h H a T -a^h. 



(32) 



(33) 



k=l 
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By noting that the pdf of h is given by fh(h) = (ir) NtNr exp(— ||/i|| 2 ), we find the estimates of r 
as follows. 

1. For finding r MAP by denoting the MAP estimate of h by h, setting x = [<r MAP ; anc l 
considering the independence of r and h, from (jlOp we get 

[f MAP ; h] = arg max f x (R \ r,h)f h (h) = arg max argmax/i(i? | r,h)f h (h), 
TeA(X),h TeA{X) h 

or equivalently 

h = argmax/!(i? | r,h)f h (h), and r MAP = arg max f\{R | r, h)fh(h). (34) 
Using (|30|) . (f33l) . and (f34l) and noting the definition of a x given in (|32|) provides that 
ft, = argmax/i(/? | r ,h)f h (h) 

h 

= argma^C^-^^+^expf- ^ ||r[/t]|| 2 + /i H A T ft - ft H a T - a^i + \\hfjj 
= argmini h H (A T + I NtNr )h - h H a T - a^h I (35) 

h v ' 



Based on the above estimate of r, D\(t, h) given in (J33]) becomes 

K 



D^r.h) = ^Hr^f-a^S- 1 ^-!^!! 2 . (36) 



fc=i 

Hence, the ML estimate of r can be found by solving 



f MAP = ar 



g max log /i(J2 | r, h)f h (h) = arg max { -Di(r, h) - \\h\\ 2 \ 



arg max {a T B T l a T \ (37) 
-reA(X) 

arg max (At + i") _1 a T } 

( N t N r , / n ,2 1 
I \ ^ \ ^Tt.nVTt-nJ 

arg max < > > v 



(38) 



which concludes the proof. 
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2. For this case which we average out the effect of h by setting x = r from (jlOp we get 

f avc = arg max fi(R I r) = arg max / f\(R \ T,h)fh(h) dh 
■reA(X) T^A{X) J h 

= arg max / exp ( — Di(r,h)) fh(h) dh (39) 

= arg max / exp I — (h H (A T + I^.^^h — h H a T — a^h) I cZAx 
^A(X)J h V V < ^— ^ V 

=-B r 

= arg ™*x J exp (- | (ft - B~ l a T ) H B T (ft - B T l a T ) - at? ' B~ l a T j) dh (40) 

= arg max exp {a^.B~ l a T ) / exp ( — { (ft — B~ l a T ) H B T (ft — Bz}a T ) \ 
-reA{X) ' J h V I v J 



dh 

= arg max exp (^a^B~ l a T ^ \Bt\~ 1 (41) 

= arg max { a 1 ? B~ l a T — \og\B T \\ 

-reA(X) 

= arg max { (A T + a T — \og\A T + I\\ 

■reA(X) 

Nt Nr ( I / \|2 /i? N 

Noting the definition of a mt n( T m,n) given in (|3ip concludes the proof. 



B Proof of Theorem [3] 

In order to find the ML estimate of r and £ form fjlQf) we have 



or equivalently 



[r, Cl = arg max f\(R \t,C) = arg max argmax fi(R I r,C), 
L J ^eA(X),( t&A{X) C 



C = argmax/i(i? I r,C) and r = arg max /x(JJ|t,t) 

C -re^(AT) 



On the other hand, from (|18|) the likelihood /i(-R | T, r) is given by 



N r K 



n=l k=l 



MR I t,C) = (tt)-**" exp ( - £ £ r„[fc] - ^/— C c"^ £ ^ r m ,, 

m=l Tm .« 



By taking into account the orthonormality assumption on s m [k;T mjn ] we get 



N r K 



n=l k=l 



D 2 (t,() = zZYl " V W ^ C X, ^ s m [A;;T m , 
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2 



(42) 



Nr Nt 



n=l m=l ' 



Nr N t rjr _p K 

J- £ e-^-M ^_ C £g_ Sm [fc;r m , r 

n=l m=l 

e i27r/ c r m ,, ; 

n=l m=l 



2? 




T TO , n k=l 



c* — 

Nt ^ J> 



r n [k] S m [k] T TO>n ] . 



(43) 



k=l 



D2(t, C) is quadratic in C and its minimum which is attained at the ML estimate of £ is given by 



( = arg max log fi(R \ T, C) = argrnin D 2 (r,() 



Z_/m=l Z—/n=l 



JZ±F **" fcTm ' n Y^UrrM s*Jk;r m , n ] 

By substituting (|44p back in (03]), and after some manipulation, D 2 (t,() becomes 



(44) 



£ 2 (r,c) = £>[fc]|| s 



E 



k=l 



T s N t 



\<\ 2 EE ( )2 • 

m=l n=l x ' ' 



Therefore, 



N t N r 

t = arg max log fi(R | r, Q = argrnin £> 2 (r, C) = arg max \ |( | 2 EE 7 \2fl f 



arg max 



3" 



' m . 77. 



(45) 



C Proof of Lemma [3] 

We provide the proof by induction. 
1. N = 2. We define 



q{tiM) = W atl gi{t 1 ) + e^g 2 (t 1 



Then 



\ti,t 2 ) = \9i(ti)\ 2 + |<7 2 (ii)| 2 + e^- t ^g 1 (t 1 )g* 2 (t 1 ) + e^~^ g^gl^), 



and therefore, 

% 2 (*i,t 2 )) 
dt 2 



ja 



g2{ti)gl{t 1 )-e^- t ^g l {t l )t 2 {t 1 ) 
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For finding max (2 q(ti,t 2 ) we set d ^ q q^ 2 ^ = and obtain 



e jat 2 = e ja tl 9l{h)g* 2 (h) 

We can check that this is a maximum by noting that 



d 2 (q 2 (h,t 2 )) 



t 2 =t2 



= (ja) 2 [e^( t "^ tl )^(t 1 ) 52 (t 1 ) + e^( tl - t " 2 ) ffl (t 1 )^(t 1 ) 
= -2(a) 2 | ff i(ti)<7 2 (ti)| <0, 



Hence, 



msxq(t 1 ,t2)=q{t 1 ,t 2 ) = Vl<7i(*i)| 2 + Mh)] 2 + 2|<7i(ti)<7 2 (*i)| = V |^(ti)|. 
t 2 < — • 

2. Inductive assumption. 



i=l 



max 

{*2,— ,tjv} 



A? 



3. Claim 



i=i 



Af+l 



JV 



= 

<9tjv + i <9ijv+i 



max 



E e ^(ti) 

i=l 

Let t = [h . . . t N+1 ] and q(t) = Ztf e jaU gi(h) 
d 



i=i 

= Ei»(*i)i- 

i=i 

, then 



JV 



i=l 



/v 



A? 



5]e^^- i -) 9 :(i 1 ) 9w+ i(fi) + ^e^( i --^) ffi (ii)ff/; + i(ii) 



i=i 



i=i 



If we set ¥^ = 0, then 



Ei=i e-^H*(ti)2AT + i(*i)' 



Also t jv+i is a maximum since 



d 2 (q 2 (t)) 



UL N+1 



tjv+i— tjv+i 



= (J«) 2 

= (J«) 2 
= -2a 2 



A 



A 



E^* +i " t ' ) A , (ii)ffw(fi) + E e3a(! '" iN+l) *( f i)^i(*i) 



1=1 



i=l 



/AT Af \ 



k=l 



N 



Y,e' aU g l {t 1 )g* N+1 {t l : 



i=i 



< 0. 
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Therefore by substituting i/v+i and some simple manipulations we get 



maxq 2 (t) = q 2 (t x , ...,t N , i N+1 ) 



N 



i=l 



+ l<?7V+l(tl)|' 



A 



N 



+ 2[Y,e> aU g t (t 1 )g* N+1 (t 1 )Y,e- jatk gUti)gN + i(h) 



vi=l 



fc=l 



A? 



i=i 

AT 

=1 

AT 



+ |(7iV + i(ti)| 2 + 2 



A 



£e**«(ti)flWfi; 



1=1 



+ |<7iV + i(ti)| 2 + 2| 5A r +1 (tO 



+ l<7iv+l(*l)l • 



A 



i=i 



i=i 



As a result, 



max q(t) 

{t2|-"!*JV+l} 



max maxg(i) = max 

{*2,."I*JV} *JV+1 {*2, •••>%} 



A 



i=i 



+ |ffAT+l(*l 



A 



AT+1 



Ei^(*i)i + i^+i(*i)i = E i«(*i> 



8=1 



1=1 



D Proof of Lemma 31 



1. By recalling the definition of S n (r) = y^'_i -, — ^—^s m e j27T ^ Tl ' 1 Tm -™) and the pdf of h given 
by fh(h) = -ie - '^ 2 , the likelihood under hypothesis Tlx is given by 



fi(R\T,h) = (vr) 



-if AV 



exp 



N r if 

EE 

n=\ k=l 



E 




r ™W ~ \lj^ S n( T )h s[k;r 1} i] 



= D 2 (T,h) 



We further simplify D^{r,h) as follows. 

AV if 



ra=l 
if 

=E 



n=l fc=l 
if 




r\k II" + \h\ x 



k=l 



E 



N r 



T s N t 



Ei^ 



n=l 



-pT JV r if r-p- AV if 

E E r « - h ]/w t S £ r » ri ^- 



71=1 fe=l 



n=l k=l 
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Therefore, the MAP estimate of h is given by 

h = argmax{/i(i? | r,h)f h (h)} = argmin \D 2 (T,h) + \h\ 2 \ 

h h <* J 



&E£i^(T)Sf =1 r w [fc]**[A:;n,i] 
ASiWr)P + l 
which in turn after some manipulations provides that 

K 



(46) 



D 2 (r,h) + \h\ 2 = ^\\r[k]\\ 2 -\h\ 2 {^-J2\S n (T 

k=l 1st n=1 



Therefore, 



.MAP 



r e Nr 

arg ^>) \ fl{R 1 T ' h)fhCh) } = arg ^') ^{tnX '^ (t)|2 + 1 

^ s n=l 



Ef=i^WEfc=i r nW s *[ fc ; r i,i] 



arg max 



which is the desired result. 



(47) 



2. By recalling the definition of S n (r) = Y^Li 7 — ^Sme^^' 71,1 Tm,n " ) and the pdf of h given 
by //j(/i) = ^e~l A l 2 , the likelihood under hypothesis Hi is given by 

/i(fl|r) = / h(R\r,h)f h (h) dh 
Jh 



exp 



N r K 

EE 

n=l k=l 



E 




r n[k] ~ \lj^ S n(r)h s[k;r lx 



+ \h\ 2 dh. 



~-D 3 (t) 



We further simplify D^(t, £) as follows. 

N r K 



^3(r,C) = EE 



n=l k=l 
K 



r n [k] - y ^-S n (r)h s[k; t 1a ] 

N r 



+ \h\' 



fe=i 



r-FT N r K rjr Nr K 

* n=l k=l 1 n=l k=l 



K 



= EHfe]ll 2 + 



k=l 



N r 



Y,\Sn(r)\ 2 + l- 



l ^ . . \/w t ^n=l S n(T)Hk=l r n[k}s*[k;T lt l} 



^T?N t Zn=l\Sn(T)\ 2 + l 
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E_ 

Nt 



Y,n=lSUT)^ =l r n lk] S *lk;T 1A ] 



T^Nt Sn.=l \Sn(T~)\ 2 + 1 



(48) 



By following the same line of argument as in the proof of Theorem [2] it can be readily 
demonstrated that the ML estimate of the time-delay for phased-array system is given by 

f ave = arg max MR I r) 



arg max < 
t<eA'{X') 



2 i Nt 



log 



E 



N r 



TJV) 



J2\S n {r)\ 2 + l 



n=l 



3. Setting P = 1 in ([3]), it is seen that h m>n = ( e J 27r /c T ™.«. Therefore, the likelihood under 
hypothesis Hi becomes 



h(R\T,() = W 



-KN r 



cxp 



/ N r K p=- N t 

-EE .w-#,.iE^ 

\ n=l fe=l m=l ' 



i2ir/c(n,l— 2r m , n ) 



By defining S^fr) = Y^'_i 7 — ^Sme* 2 "^^ 1 ' 1 2r m , n ) anc j f n wing the same line of argu- 
merit as in the previous part the ML estimate for £ is found as 



C = argmax{log/i( J R | r, C)} = arg mm {£> 4 (t,C)} 



and consequently 



A 



d 4 (t,c) = £>[*]|| s 



E 



T s N t 



ici 2 Ei^ 



k=l " u n=l 

Therefore by following a similar approach as in part 1 we get 



f E 

f = arg max fi (R I r, C) = arg max < — 



Nr 



= arg max 
which is the desired result. 



££i3t(T)££ =1 r n [fc]s*[fc;n,i 



A" 



n=l 
|2 



(49) 



E Proof of Lemma [6] 

We define the events 

M r M ~\ M 

Vi = |J {l£ < 7} , and P 2 = J]^< 7 , and P 3 = |J {l£ < ^} 

m=l Lm=l J m=l 
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Clearly as P 3 CD 2 C V\ we have Pr(£> 3 ) < Pr(£> 2 ) < Pr(£>i), where Pr(£> 2 ) is the probability of 
interest in (|2ip . In order to analyze P(T>i) we define random variables 7 m ~ AA(p • |/i m |,cr 2 ), for 
which we have 

Prftf < 7) = PrC«<7HPr(-v^<^< 

= q( ^2 ) (50) 

= q(^), (51) 

where the transition in (|50p holds by noting that for any choice of 7 > a 2 , from the table of Q(-) 
function it is readily verified that for sufficiently large values of x, and any fixed value A > 0, 
Q{x) 3> Q(x — A). The asymptotic equality in (f5"Tj) is justified by noting that p\fJ, m \ 3> yfi. By 
further using the following known bounds on the Q(-) 



Vx G M + , — L- e~ x2/2 < Q(x) < 

V2^x\ x 2 J ^x 



v 27rx \/2ttx \/2irx 



we get 



By further taking into account that J- e x2 / 2 = e a;2 ' /2 and noting Remark Q] we get 

Q(x) = e - 2 / 2 . (52) 

Invoking the exponential equivalence in (|52p on (|5ip provides 

2 2 

Pr(K^< 7 ) = e-^. (53) 
As /i m is an unknown random variable distributed as AA(0, cr 2 ), by averaging out its effect we get 

E Mm [Pr(Y^ < 7 )] = — / e ~ ( ^ + ^ - ^ m = _ = p ~\ (54) 



Next, by considering the statistical independence of {Y m }M =1 we consequently find 

M 

E M [Pr(Px)] = [] E Mm [Pr(y, 2 < 7)] = P" M (55) 

m=l 

Note that, as long as the asymptotic behavior is concerned, the asymptotic exponential order 
of E M [Pr(2?i)] does not depend on the choice of 7 and by following the same line of argument as 
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above we can extend the same result for the choice of 44 and thereof for the exponential order of 
E M [Pr(X>3)]. As a result, we can similarly show that 

E M [Pr(P 3 )]=/^ M (56) 
By noting that Pr(2?3) < Pt(T>2) < Pr(£>i) and using ((351) and (156]) we get that 

Pr(P 3 ) < Pr(^ 2 ) < Pr(Pi) ^ Pr < J\ = Pr(V 2 ) = p~ M , 

\m=l / 

which concludes the proof. 
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